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ABSTRACT 

At least 50% of the baryons in the local universe are undetected and predicted 
to be in a hot dilute phase (10 5 -10 7 K) in low and moderate overdensity environ- 
ments. We searched for the predicted diffuse faint emission through shadowing 
observations whereby cool foreground gas absorbs more distant diffuse emission. 
Observations were obtained with Chandra and XMM-Newton. Using the cold 
gas in two galaxies, NGC 891 and NGC 5907, shadows were not detected and 
a newer observation of NGC 891 fails to confirm a previously reported X-ray 
shadow. Our upper limits lie above model predictions. For Local Group studies, 
we used a cloud in the Magellanic Stream and a compact high velocity cloud 
to search for a shadow. Instead of a shadow, the X-ray emission was brighter 
towards the Magellanic Stream cloud and there is a less significant brightness 
enhancement toward the other cloud also. The brightness enhancement toward 
the Magellanic Stream cloud is probably due to an interaction with a hot ambient 
medium that surrounds the Milky Way. We suggest that this interaction drives 
a shock into the cloud, heating the gas to X-ray emitting temperatures. 

Subject headings: Galaxy: halo — intergalactic medium — Local Group — Mag- 
ellanic Clouds — galaxies: individual (NGC 891, NGC 5907) — X-rays: diffuse 
background — X-rays: galaxies 
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Introduction 



An inventory of the baryons at the present epoch shows that the visible galaxies and 
the hot gas in clusters and rich groups compri se only a fraction o f the baryons inferred 
from cosmic microwave background observations ( Spergel et al. 12007 ) or fro m the D/H ratio, 
combined with Big Bang nucleosynthesis calculations (IKirkman et al.ll2003l ); these two values 
are in good agreement. Most of the baryons must be in a form that is difficult to detect, 
and cosmological simulations suggest that it is still in the gaseous phase in regions of low 
and moderate oyerdensities (p/ A p = 10 — 100; iFukugita et al.l Il998l ; ICen fe Ostrikerl Il999 



Dave et al.ll200ll : ICroft et al.l 1200 ll ; ICen fe Ostrikel 120061 ) . In the simulations, gravitational 
collapse leads to compressional and shock heating of the gas, producing temperatures of 10 5 
- 10 7 K, where hydrogen is almost completely ionized. The metallicity of this gas increases 
with cosmic time due to enrich ment by galactic winds, leading to present day metallicities 
of about 1/10 of the solar value (IDanforth fe Shulll 120051 ) . 



At X-ray energies, one can try to detect this gas in emission or in absorption and 
each approach poses substantial challenges. For absorption line studies, the intergalactic 
medium (IGM) is predicted to produce only small equivalent width lines that are at the 
limit of detectability (reviews by iBregmanl 120071 ; iNicastro et al.ll2008l ). Another approach is 
to search for emission from the IGM, since the presence of metals, and especially of Fe and O, 
produces X-ray line emission, which when summed over many filaments and large volumes , 
should be detectable as a diffuse X-ray background (ICen fe Ostrikerlll999l ; ICroft et al.ll200ll ). 
This diffuse background is near the present limit of detectability and is predicted to be less 
than the typical instrumental backgrounds or than the point source background, even in the 
0.4-1 keV range where IGM emission is expected to be most prominent. It is not possible 
to subtract the unresolved point source background and the instrumental background to 
sufficient accuracy, as they can vary from field to field and are not known with sufficient 
accuracy a priori. 

A detection method that avoids this problem is that of X-ray shadowing, whereby an 
excess of cold gas (e.g., a gas cloud) lies in front of the emitting source and falls within the field 
of view of the detector. The cloud absorbs some of the background X-ray emission, producing 
a local minimum on the image. To apply this method, it is crucial that the distance to the 
absorbing cloud is known and that it's absorption column has been measured independently, 
such as from H I measurements. This shadowing technique has been app lied within the 
Milky Way to show that there is a hot halo above the disk of the Galaxy (ISnowden et al. 



19971 ) and that there is X-ray emission from the Galactic bulge (IPark et al.lll997l ). 



Here we determine whether shadows are cast in two situations that address the missing 
baryon problem. The first set of observations uses the H I in external galaxies to shadow 
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the bac kground emission from th e ensemble of hot cosmic filaments. This was utilized in the 
work of iBregman fc Irwin! (120021 ). where the dense atomic and molecular gas in the disk of 
NGC 891 appeared to produce a shadow at the 98% confidence level. That is, in the 0.4- 
1.0 keV band, there were 39.5 photons in the extinction region compared to 57.6 photons 
in the non-extinction region, a deficit of 18 photons. Since this is less than a 3a effect, 
it is important to confirm or refute the result, and as the initial observation was 51 ksec, 
there is room for improvement. A significantly longer observation (120 ksec) was obtained 
for another program, which permits us to examine the quality of the shadow seen with the 
shorter observation. A second shadowing test is also presented, whereby we use the very 
extended H I from the edge-on galaxy NGC 5907 to search for a shadow against the diffuse 
IGM emission. 

The other location that one might hope to find the missing baryons is in the Lo- 
cal Group, where the X-ray absorption line s from O VII and O VIII have be en detected 
( jNicastro et al.ll2003l ; iRasmussen et al.l 120031 ) as has OVI (ISembach et al.l 120031 ) . These ab- 
sorption lines are at zero redshift, so the gas could either be a Galactic halo or a Local Group 
medium. If it is a Galactic halo (typically 50 kpc or smaller), the mass involved i s less than 
t he known mass of cold gas in galaxies, and this is argued for by IWang et al.l (120051 ) and 



Bregman &: Lloyd-Daviesl (120071 ). However, if t he gas fills the Local Group, then the hot gas 



mass can dominate the total baryon content (jNicastro et al.l 120021 . 12003 



Rasmussen et al. 



20031 ). as expected from the missing baryons. Since t here is also diffuse X-ray emission that 
is associated with this gas (IMcCammon et al.ll2002l ). one can hope to use the shadowing 
technique, provided that there are absorbing clouds at distances of at least tens of kpc from 
the Milky Way. There are such clouds, as the Magellanic Stream has some dense knots 
of H I and it lies at about 50-60 kpc. Anoth er class of potential targets are the compact 
high velocity clouds, which iBlitz et al.l (119991 ) suggested to be orbiting through the Local 
Group at distances of ~ 10 2 5 kpc. However, new evidence suggests that these are associ- 



ated with the Galactic halo and lie within 50 kpc (IMaloney fc Putmanll2003l ; iPutman et al. 



20031 ; IWestmeier et al.l 120071 ) . We have obtained observations toward these clouds suitable 
for shadowing studies and report on them below. 



2. Observations 

There are five different observations discussed below of two edge-on galaxies and two 
H I clouds in the Local Group. As some observations were obtained with Chandra while 
others obtained with XMM, we will discuss the results as a function of the instrument used. 
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2.1. Chandra Data and Data Reduction 

The Chandra observations that we analyze in this work were obtained with the ACIS 
instrument with the objects centered either near the aimpoint of the S3 chip or near the 
center of the S3 chip. Basic observational information on the individual pointings are given 
in Table 1. 

For the data reduction, which has become fairly standard, we used CIAO 3.2.2 along 
with CalDB 3.1.0. In each set of data, we identified the bad pixels and pixels affected by 
cosmic ray afterglows and excluded them from the data set; only grades 0, 2, 4, and 6 were 
used. Periods of high background were excluded, with the threshold being 40% above the 
modal value of the background. Exposure maps were produced for the S3 chip, on which 
our analysis was focused, and the images were either flattened or the exposure map was 
used separately to correct the extracted counts. For the identification of point sources, we 
used both wavdetect and direct identification by eye, which led to nearly identical results. 
The observations of the Magellanic Stream Cloud were obtained in FAINT mode while the 
observation of NGC 891 and the Compact High Velocity Cloud were obtained in VFAINT 
mode. 



2.1.1. Magellanic Stream Cloud MS30. 7-81. 4-118 



One of the best targets is a neutral gas cloud at a distance near the edge of what might be 
considered the Galactic Halo, about 50-70 kpc. At this distance is a set of clouds, the Mag- 
ellanic Stream (IGardiner fc Noguchilll996l ). which passes over the South Gal actic Pole. Re 



cently, a deep survey was obtained with the Parkes 21 cm multibeam system (IPutman et al. 



20031 ). revealing a great deal of structure in the Stream at 15' resolut i on (Fi gure ED). Higher 
resolution (1') radio synthesis observations were obtained by iBriind (119981 ) of the highest 
column density regions and one region shows structure with H I columns exceeding 4xl0 20 
cm -2 (Figure [2]), significantly greater than the Galactic H I column of 1.5x10 



20 cm 2 . 



Two partly-overlapping images were obtained, on the H I maximum, and away from the 
H I maximum. The exposure on the H I maximum was one of the rare times when there were 
no flares and the background was relatively low, so nearly the entire exposure was used and 
this was the most useful of the two exposures. Potential shadowing by the Magellanic Stream 
cloud should have the most pronounced effect at energies below 1 keV, so we produced a 
flattened 0.4-1.0 keV image with point sources removed (they accounted for 20% of the total 
counts in the image). In this image, there is an apparent extended symmetrical brightening 
at 00:13:32.6, -27:11:07 (RA, Dec, J2000), which may be a background group or cluster; it 
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appears in both the soft and hard (1-6 keV) bands. 

To search for shadows, four regions were chosen, from the highest to the lowest H I 
regions. Using the flattened hard image, all four regions were consistent with a constant 
value for the surface brightness, being within 1.3cr of the mean. In the 0.4-1.0 keV image, 
there is some modest variation, with the high H I column regions being brighter than the 
low H I regions, the opposite of the shadowing signature. When we compare the two highest 
H I regions to the two lowest H I regions, the surface brightness values are 1.29 ± 0.04 
cts ksec -1 arcmin -2 and 1.20 ± 0.04 cts ksec -1 arcmin -2 , a difference of 0.093 ± 0.058 cts 
ksec -1 arcmin -2 . This difference corresponds to a modest brightening of 1.6a, rather than 
the dimming that a shadow would produce. 

There will be fluctuations in the X-ray sky due to variation in the Galactic X-ray 
background, which we can assess in this part of the sky by using our off-cloud data. Even 
with a longer observation, it may not be possible to improve the limit on the shadow due 
to the fluctuations in the soft Galactic X-ray background (this reg ion has a fairly smooth 



background, as seen from the ROSAT data (ISnowden et al.lll997l ). We analyzed the off- 



cloud images, which after screening high background regions (more than 50% higher than 
the modal value), led to a useful observing time of 33.49 ksec. After flattening the image, 
we divided it into nine regions about 300 pixel square, about 15% smaller than most of the 
regions used for the on-cloud study. The errors per bin were about 6% and the 0.4-1.0 keV 
background in the field is consistent with a constant value, as the variation is ±1.4cr. 

Next, we need to separate the amount of flux from the instrument and from the X-ray 
background. The instrumental background for the ACIS S3 chip was fairly stable over the 
past five years (since 2000.5) and lower than at launch. The total chip background for the 
observations is expected to be fairly modest (and similar) for the observation dates of our 
targets (7.5-9 cts s -1 chip -1 ; see Figure 6.24 in MPOG v8). Using the standard grades, 
avoiding flares, and for the 0.4-1.0 keV band, the instrumental count rate is expected to be 
about 4.2 xlO -8 cts pixel -1 for 0.5" pixels, or 0.61 cts ksec -1 arcmin -2 . This corresponds to 
about half of the counts registered in the above measurements. 

Another component that needs to be subtracted is the unresolved point source contri- 
bution from AGNs. For a 50 ksec Chandra ACIS-S observation, where the limiting point 
source is defined by five photons, the unresolved point source contribution to the diffuse 
emissi on is about 15% of the point source component of the background in the 0.5-2.0 keV 



band (IMoretti et al.ll2003l ). When converted to the 0.4-1.0 keV band, the unresolved point 
source contribution is 0.037 cts ksec -1 arcmin -2 . This unresolved point source contribution 
is less than 10% of the diffuse background for all Chandra observations and closer to 5% for 
most of the Chandra observations. We remove this small component, correctly scaled to the 
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observing time, in our analysis. 

We remove these various components to obtain the diffuse emission from gas in the 
shadow and non-shadow region. With minor differences, this procedure is followed for the 
other observations. There is one final consideration, which is the removal of the emission from 
the Local Hot Bubble (LHB). Once this is removed (discussed below), we have measurements 
of true diffuse emission from above the disk. 



2.1.2. CHVC 125+41-207 

Most of the high velocity clouds are resolved easily with single dish radio telescopes, 
but there are a number that have significant columns yet have s mall angula r sizes. These are 



the compact high velocity clouds (CHVC). iBraun fc Burton! (119991 . |2000| ) mapped several 



CHVCs and a few of them have column densities of up to 4xl0 20 cm -2 , with sizes of 10-50 
arcmin 2 . These clouds are ideal targets for X-ray shadowing since they are mapped at good 
S/N in H I and have high density regions that fill a significant fraction of a single ACIS chip. 

In selecting a CHVC as a target, we needed a cloud whose column density is comparable 
to that of adjacent regions in the Galaxy and whose solid angle covers a significant fraction 
of the ACIS S3 chip, which is the most advantageous for this prog ram. Of the clouds tha t 



have been mapped with radio synthesis arrays in the 21 cm H I line (IBraun fc Burtonll2000l ). 
the cloud CHVC 125+41-207 is the most suitable (see Figure EJ). It is at I, b = 125°, 41° 
(and velocity -207 km s _1 ), with the highest column density part confined to an elongated 
structure of dimensions 6' x 17'. At lower column density levels, si ngle-dish data show that 



it has a head-tail configuration (Figure H]) and is about 0.8° by 2.5° ( iBriins et al.ll200ll ). The 
synthesis H I map, which has 28" resolution, shows that the highest column density part 
of the cloud is about 5xl0 20 cm -2 , which is larger than that of smooth regions adjacent to 
the cloud (3xl0 20 cm -2 ), which are due to the Galactic disk. Aside from the usual point 
sources, there is a previously unknown diffuse feature that occurs just to the east of the 
H I cloud with an RA and DEC (J2000) of about 12:29:38.7, 75:21:00, and with a radius of 
about 30". This falls in the background subtraction region of the image and was excluded 
when quantifying the shadow. There is no obvious optical counterpart on the Digital Sky 
Survey images, so it is most likely either a distant cluster or it is caused by the interaction 
of the CHVC with its surroundings. 

This cloud was observed three times at different positions and roll angles. In all cases, 
we identified and excluded point sources from the analysis and defined on-cloud and off- 
cloud region. In the first observation (Obsid 2484), there were several long periods of high 
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background, so the final useful data was 10.1 ksec. In the second observation (2486), the 
instrumental background was well-behaved and we were able to use most of the data, 21.2 
ksec. Conditions were most favorable in the last observation (2253), which is twice as long as 
the preceding observation and the orientation of the chip is ideal so that about half contains 
high-column H I and half has the lowest H I column of any chip position; the useful time 
was 45.3 ksec. The three observations yield a useful integration time of 76.5 ksec, so we 
differenced the regions of higher H I column (5xl0 20 cm -2 ) and lower H I column (3xl0 20 
cm -2 ). The mean surface brightness in the higher H I regions is 1.20 cts ksec -1 arcmin -2 , 
which was 4.9% brighter than the surface brightness of the lower column density region. As 
the uncertainty in the difference is 3.8%, the brightening is only a 1.3<r effect, although in 
the opposite sense of the shadow that we were searching for. 



2.1.3. NGC 891 

The more recent and longer observation of NGC 891 had only a few periods of high 
background so that 108.5 ksec were used out of the total exposure time of 120.4 ksec. Not 
only is this useful time a factor of three greater than the 36 ksec used in the earlier exposure, 
but the instrumental background appears to be lower. In addition, two strong sources in 
NGC 891 are seen in both exposures (one source is SN 1986 J), so a very high relative accuracy 
is achieved and we can be certain that the regions used to analyze the earlier observation 
are sampling the same J2000 coordinate regions. 



The shadow region was defined in lBregman fc Irwin! (120021 ) by the dark extinction region 



that obscures nearly all of the optical light from parts of the disk. All point sources were 
excluded, an important consideration because there are several point sources detected from 
the disk, most of which are probably X-ray binaries. The extinction region, with the point 
sources masked out, constitutes 1378 arcsec 2 (see Fig. [5]). For the background region, we 
needed to obtain a surface brightness that would not be contaminated significantly by the 
extended emission from NGC 891, so these regions are in the outer part of the S3 chip, 
perpendicular to the disk of the galaxy. The exposure corrected background on the SE 
side is 1.22 ± 0.04 cts ksec -1 arcmin -2 and on the NW side it is 1.39 ± 0.04 cts ksec -1 
arcmin -2 , so there may be a slight gradient in the Galactic background; the average of the 
two backgrounds is 1.31 ± 0.08 cts ksec -1 arcmin -2 . In comparison, the surface brightness 
in the extinction region is 1.80 ± 0.24 cts ksec -1 arcmin -2 , a difference of 0.49 ± 0.25 cts 
ksec -1 arcmin -2 , which is 2a above the background level, whereas the earlier observation 
had found it to be about 2cr below the background level. This slight excess is also present 
when we subdivide the energy bands into 0.4-0.7 keV and 0.7-1.1 keV (about a 1.4a excess 
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in both bands). Thus, the shadow previously seen in NGC 891 is not confirmed and we 
conclude that the earlier observation was a statistical fluke. 

2.2. XMM Data and Data Reduction 

2.2.1. NGC 5907 

We have examined the H I maps for all edge-on galaxies as well as the gas torn away in 
interacting galaxies and we found that the best combination of optical depth and solid angle 
is given by the H I in NGC 5907. The gaseous disk of NGC 5907 not only extends beyond 
the stellar disk, it has a modest warp where the gas turns away from the smaller stellar disk 
(Figure E]). The column density of this extended gas rises above 1 x 10 21 cm" 2 , which causes 
significant absorption in the 0.4-1.0 keV region. 

The observations of NGC 5907 were taken February 20 and 28, 2003 (observations 2 
and 1, respectively, with total exposure times of 52.74 ksec and 52.19 ksec). The data 
were extracted using the Science Analysis Software (SAS) version 7.0.0 of the XMM-Newton 
group. Unfortunately, the observations were contaminated by strong flaring. We excluded 
the time intervals of the flares and in addition applied the standard three-sigma clipping to 
the data based on the 12-14 keV light curves. We checked the resulting good time intervals 
in the lower energy band (0.4-7 keV). No additional time intervals had to be excluded. Due 
to the strong flaring, only 12.5 ksec of exposure time were usable for NGC 5907 (about 12% 
of the total exposure time). For our analysis, we used only the PN data because of the higher 
sensitivity of that CCD. 

The data were binned by 50 pixels during SAS processing. The lowest energy range 
(0.2-0.4 keV) is dominated by an instrumental background, whereas the background in the 
highest energy range (1.0-10.0 keV) is expected to consist mostly of unresolved AGNs. We 
therefore extracted a 0.4-1.0 keV energy band image for our shadowing experiment. In 
addition, a 0.5-2.0 keV and a large energy band image (0.2-10.0 keV) were also created to 
aid in our analysis. 

In order to identify point sources we created two smoothed images for either observation 
by applying Gaussian filters of 5 and 25 binned pixels, respectively, to the 0.2-10.0 keV 
images. The difference images between these two smoothed images revealed point sources 
and small scale structures from which we created a bad pixel mask for either observation. 
Remaining possible point sources were marked by hand and included in the bad pixel mask. 
Dead columns were also identified and masked. CCD boundaries and hot pixels were derived 
from the exposure maps and included in the bad pixel masks as well, which were then applied 
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to the filtered images. 

In this shadowing experiment, a flat background is crucial to obtain reliable measure- 
ments. While part of the X-ray background is caused by scattered high energy particles 
creating a smooth background Bd at , high energy particles entering the aperture and follow- 
ing the optical path are subject to the vignetting effects of the mirrors creating a focused 
background B {oc . We attempted to disentangle these two components in each energy band 
by measuring the background counts G\ (i = 1-15) in 15 circular aperture (each about 100" 
in diameter) spread across the NGC5907 observation away from a possible galaxy X-ray 
halo. We also obtained the exposure times t\ for these regions using the exposure map. The 
flat background then could be derived for each set using the least \ 2 method on 

Ci = -Bfoc ■ ti + -Bflat- 



We subtracted the B& at component, divided by the exposure map and normalized the 
image by multiplying it with the actual exposure time of 12.5 ksec. The resulting flattened 
image displays an apparent inverse vignetting pattern (e.g., positive pixels values increase 
towards the outer parts of the field of view as shown in Fig. [7]). To understand this pattern, 
one has to consider the fact that most pixels of the X-ray background in our data did not 
receive any photon events even after the binning. The vignetting pattern in the data is 
therefore not represented by an increase in pixel values toward the optical axis, but by an 
increase in photon event density over a larger area. In the flattened image, the former 
non-zero pixel values therefore increase toward the edges of the CCD, whereas the former 
zero-value pixels are decreasing to more negative values at the edge of the field of view. 
The larger density of negative pixels thus compensates the increase in positive pixel values 
toward the edges of the flattened image. 

For extended sources filling the field of view (i.e., MS30.7), this method of deriving the 
flat background component is not possible. The basic vignetting pattern cannot be removed 
from observations of extended sources (this is a more significant problem than in the Chandra 
observations that used only a single chip with a smaller field of view). W e ther efore tried 



a second method using the closed filter wheel (CFW) data of iMarty et al.l (120031 ) to derive 
the flat background component. For this CFW method, we chose measurement regions 
along contours of constant effective exposure time to exclude the vignetting effects from 
our measurements. The instrumental background shows sp atial inhomoKeneities at certain 



energies due to the electronics board behind the CCD (e.g., iFreyberg et al.l I2004J ) . Most of 
these fluorescence lines, however, are at higher energies and thus should not affect the 0.4-1.0 
keV band. We measured the instrumental background in the same regions of the CCD as 
for NGC 5907 and subtracted those count rates from our galaxy measurements accordingly. 
For better comparison, we chose the same measurement regions for both methods. 
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We used the H I contours of IShang et al.1 (119981 . our Fig. [6]) to determine the area 



that could absorb X-ray background photons. The area of the galactic disk has the highest 
absorbing column, but it is also more likely to be contaminated by unresolved X-ray point 
sources. We therefore chose regions in the H I halo of possible shadowing and regions 

outside the halo as X-ray background. We tried to make the measurement regions as large 
as possible while minimizing the number of masked pixels in each region (as those may have 
different effects on the count rate due to the vignetting). Figure [8] shows the two region 
pairs, one in the northern and one in the southern part of the field of view. The number of 
actual pixels in each region was obtained from the bad pixel mask to correct the measured 
photon counts. The measurements of both regions are listed in Table [2] for both methods. 

Regions 1 and 2 yielded on-minus-off-galaxy measurements of —0.45 ±0.21 and +0.29 ± 
0.23 cts ksec -1 arcmin -2 , respectively, for the x 2 method and —0.60 ±0.34 and ±0.34 ±0.39 
cts ksec -1 arcmin -2 for the CFW method. This indicates a deficit in the sense of a shadow 
at the 2cr level in the southeastern part of the galaxy halo, whereas the northwestern halo 
shows no shadow. Assuming a smooth X-ray background in the northwestern halo (and thus 
averaging the on- and off-measurements), we derive a 3a upper limit of 21-26% for a shadow 
in this region depending on the method used for the flattening. 

The detection limit for point sources in our combined observations of NGC5907 in 
the 0.5-2.0 keV ba nd is 3.6 x 10 -15 ergs s -1 cm -2 . Based on the source distribution of 



Moretti et al.1 (120031 ) in the 0.5-2.0 keV band and our detection limit, we derive a contribution 
of unresolved AGNs and galaxies of about 33% in our observations. This corresponds to a 
background fluxS for unresolved sources of about 0.25 cts ksec -1 arcmin -2 in the 0.4-1.0 
keV band. When this is removed from the potential shadow, it becomes a statistically 
insignificant la fluctuation. 



2.2.2. MS30. 7-81. 4-118 

The observation of the Magellanic Stream cloud MS30.7 was taken January 3, 2004. 
The data extraction and calibration was the same as for NGC5907 described above. The 
light curves of MS30.7 showed little flaring, resulting in 37.0 ksec of good time intervals 
(71% of the total exposure time). Due to the extent of the H I cloud on the XMM field of 
view, we did not try to derive the smooth background component using the least \ 2 method 
described above. We again chose measurement regions that followed constant exposure time 



1 Conversion between count rates and fluxes was performed in WebPIMMS at 
http: / /heasarc.gsfc.nasa.gov/Tools/w3pimms.html| using a power law index of 1.4. 
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contours of the exposure map to exclude vignetting effects and applied the CFW method 
to measure the background. Our first measurements yielded an X-ray excess on part of the 
cloud. We therefore mapped the X-ray emission along the segments of a circle centered on 
the exposure map peak (Fig. [9]). Table [3] lists the individual fluxes of all 20 regions and 
the flux differences relative to region 6, which has the lowest flux of all the regions. Fig. [10] 
displays these flux differences to emphasize the X-ray excess on the eastern part of the cloud. 
The area where the brightening occurs is consistent with that found in the Chandra data. 
Using regions 1-3 for the X-ray bright area and 6-11 as off-cloud, low H I area, we derive 
average fluxes of 2.54 ±0.09 and 1.90 ±0.05 cts ksec -1 arcmin -2 or a difference of 0.64 ±0.10 
cts ksec -1 arcmin -2 . This corresponds to a 6.4cr brightening, which is larger than the 1.6cx 
increase measured in the Chandra data. The main reason for the difference between the two 
data sets is the field of view with the XMM-PN data (especially regions 6-11) reaching the 
edge of the cloud, whereas the Chandra field of view is situated fully on MS30.7. 

The detection limit for point sources in the 0.5-2.0 keV b and of our observatio n is 



2.25 x 10 -15 ergs s -1 cm -2 . Using the distribution of sources by iMoretti et al.l (120031 ). we 
derive a fraction of unresolved sources for the background of 28%, corresponding to 0.21 cts 
ksec -1 arcmin -2 in the 0.4-1.0 keV band. The X-ray brightening of the Magellanic Stream 
cloud could therefore be as high as 0.85 ± 0.10 cts ksec -1 arcmin -2 . The 3a upper limit for 
a shadow based on regions 6-11 corresponds to 8% of the off-cloud background, which is 
consistent with the range derived from the Chandra data. 



3. Interpretation of the Results 

3.1. The Absence of X-Ray Shadows 

Our observations failed to show a shadow cast by the higher column density regions, 
although the observation for NGC 5907 produced an order of magnitude less data than 
expected, so this observation does not provide strong constraints. For our two best cases, 
MS30.7 and CHVC 125, the pointings toward the higher H I column regions showed a slight 
brightening rather than a shadow (Table 4). For CHVC 125, this brightening is only at the 
1.6er level, and the uncertainty in the ratio of the two pointings (on-cloud and off-cloud) is 
8.5%, or a 3a upper limit to the shadow of 25%. For MS30.7, the observations with Chandra 
and XMM-Newton are consistent in that the flux differences are the values expected from 
the difference in instrumental sensitivities. The difference in the ratios of the on-cloud and 
off-cloud fluxes is about at the 2a level, and if not merely due to statistical fluctuations, this 
may be caused by the slightly different regions used by the two instruments. The field of 
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view is larger for the EPIC camera on XMM-Newton than for the ACIS S3 chip on Chandra, 
so it was possible to take off-cloud regions further from on-cloud regions, which may have led 
to the greater contrast. In these observations, the flux is brighter for the on-cloud pointings, 
with a S/N of 1.7 for the Chandra measurements and 7 for the XMM-Newton data. When 
combined by weighting, the brightening is about a 7a feature. This could occur if there is 
an interaction between this Magellanic Stream cloud and a hot diffuse medium. Ignoring 
the brightening, the 3a detection limit is 14% of the diffuse background. 

In forming these ratios and upper limits from the on and off-cloud pointings, we have not 
removed the contribution from the Local Hot Bubble, which is certainly not being shadowed 
by the objects that we pointed at. We will assume that the contribution from the Local Hot 
Bubble is the same in the on and off-cloud locations, in which case, that contribution must 
be estimated. 

To estimate this contribution from the Local Hot Bubble in the 0.1-0.28 keV region 
(ROSAT bands 1+2) we assume it is equal to th at of the contribution above the disk, which 



is similar to values discussed in the literature (ISnowden et al.ll2000l ). There is significant 
variation in this ratio around the sky and there are regions where the halo contribution 
dominates that from the LHB and regions where the reverse is true. The field for the CHVC 
(I, b = 125°, 41°) is in a fairly typical region for the northern Galactic hemisphere, where the 
1/4 keV background is intermediate in intensity between the bright regions near the poles 
(the ROSAT measurements of the background at this location are R12 = 638 and R45 = 137 
in units of 10~ 6 cts s _1 arcmin -2 ). Separation between the Local Hot Bubble and emission 



above the disk can be made only on statistical grounds in general and ISnowden et al.l (120001 ) 
show that in the Northern Hemisphere, the two components are comparable. For the field of 
MS30.7, which is very close to the South Galactic Pole (I, b = 30°, -82°), the background is 
very similar to the CHVC 125 (R12 = 605 and R45 = 128 in units of 10" 6 cts s" 1 arcmin" 2 ), 
so we assume the same ratio of emission from the LHB to that above the plane. 

The 1/4 keV emission from above the disk is absorbed by the disk gas, but this absorp- 
tion decreases with energy so that in the 0.4-1.0 keV band, the emission from above the 
disk becomes relatively greater. Assuming that the temperature and abundances of the two 
components are the same (about 0.1 keV), 70-75% of the diffuse 0.4-1.0 keV background 
originates above the disk. If the component above the disk is hotter as well, the fraction is 
greater. When we assume that 25% of the diffuse emission in the 0.4-1.0 keV band is due to 
the Local Hot Bubble, there is an increased fractional error (by 33%) in the remaining flux, 
so the 3a limits for MS30.7 and CHVC 125 rise to 18% and 34% (Table gj). 

For CHVC 125, the mean difference between the higher and lower H I column density 
regions is 5xl0 20 cm -2 versus 3xl0 20 cm -2 . If most of the emission above the disk lay 
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beyond the CHVC and its gas temperature is 0.1 keV, the shadow would produce a deficit 
of 17.1% (for a temperature of 0.2 keV, that reduces to 14.3%). This would have been 
about a 1.5-2er shadow rather than the observed 1.7a brightening. The failure to see a 
shadow for this object does not place strong constraints on the nature of the Galactic Halo 
or Local Group med i um. Since this CHVC probably lies within 1 00 kpc of the Milky Way 
( jPutman et al.ll2003l ; iMaloney fc Putmanll2003l ; iPisano et al.l 12004 ). it could interact with a 
hot Galactic Halo to produce X-ray emission that could exceed the deficit expected from 
a shadow. In addition to the marginal brightening of X-rays toward the cloud, there is 
an X-ray enhancement just east of the cloud edge that may signal an interaction of the 
neutral cloud with a hot ambient medium. This is not the first time that a modest X-ray 
brighten ing has been seen toward or near a HVC, as such association s have been reported 
upon by iHerbstmeier et ajj (119951 ) , iKerp et al.l (119961 ) , and iKerp et al.l (Il999l ) . 



The shadowing constraints are a bit stronger for MS30.7, where the fractional uncer- 
tainty is 6%. The difference between the higher and lower column density H I regions is 
4xl0 20 cm -2 compared to 2.5 xlO 20 cm -2 , and for emission from a 0.1 keV thermal plasma, 
the shadow depth would be 13.1%, or a predicted 2.2cx shadow. In contrast, we find a bright- 
ening above the 7a level, so the emission from the cloud obscures the predicted shadow, 
should one exist. The most likely explanation for the brightening is that this cloud, at a 
distance of 50-70 kpc, is interacting with a hot ambient environment (see below). 

The failure to detect a shadow of the diffuse cosmic X-ray background in the Chandra 
data does not place a very useful upper limit on this component. The limit is that less than 
50% of the diffuse X-ray background comes from beyond MS30.7 and CHVC 125. However, 
if one could detect excess X-ray emission from several clouds with known distances around 
the Milky Way, it might be possible to probe the density of this hot Galactic halo. 

The observation toward NGC 5907 was designed to determine the fraction of the X-ray 
background that could be due to a cosmological diffuse component arising from the warm-hot 
intergalactic medium (WHIM). The surface brightness of the WHIM is model-dependent, 
but it could be as la r ge as 20% of the soft X-ray background in the 0.4-1.0 keV range 
( ICen &: Ostrikerlll999l ; ICroft et al.ll200ll ). Unfortunately, the measurement through the H I 
of NGC 5907 only yielded a la rms of 17%, so this observation does not provide an important 
constraint. If the X-ray background from the WHIM is near 20%, then an observation that 
is an order of magnitude longer (the original observing request), the rms would drop to 6% 
and a shadow would be seen at the 3a level. However, if the emission from the WHIM is 
closer to 5% of the soft X-ray background, the time required to detect the shadow would 
exceed 1 Msec with current instruments. 
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3.2. The X-Ray Enhancement Around MS30. 7-81.4-118 

One of the positive outcomes of this work is the detection of an X-ray enhancement at 
the location of the Magellanic Stream cloud MS30. 7-81. 4-118. From the measurement with 
EPIC, the net excess flux toward the cloud is 9xl0~ 15 erg cm -2 s -1 , for an assumed thermal 
spectrum with kT = 0.3 keV and 20% solar metallicity (the flux is not strongly dependent 
on the metallicity). For a distance to the cloud of 60 kpc, the implied luminosity is 4xl0 33 
erg s~ x . Also, the effective diameter of the cloud is about 70 pc, but it is not spherically 
symmetric. 

This cloud presumably sits in the ambient medium that surrounds the Milky Way, but 
whose parameters are poorly known. The density of the ambient medium 60 kpc from the 
Milky Way is likely to be near ~ 10 -4 cm -3 , and it cannot be larger than a few times this 
value or it woul d exceed some of the d ispersion measures toward Large Magellanic Cloud 
(LMC) pulsars ( [Manchester et al.ll2006l ). Also, significantly larger electron densities along 
60 kpc path lengths are rule d out by both the OVH absorp tion line measurement and the 
diffuse soft X-ray emission (IBregman fc Lloyd-Davied 120071 ) . A lower limit to the mean 
ambient density is given from an argument i nvolving the ram-pressur e stripping of materia l 
within dwarf sphe r oidals in the Local Group (IBlitz fc Robishawil2000l ; iGrcevich et al.ll2009l ). 
Blitz fc Robishawl (120001 ) argue that within 250 kpc, the mean gas density must be at least 



2.5 xl0~ 5 cm -3 . The temperature of this gas is likely to be about 3xl0 6 K, if it is near virial 
equilibrium with the potential well of the Milky Way and the Local Group. This would 
imply that the ambient pressure is P/k ~ 10 2 ' 5 (n/10 -4 c?r?r 3 ) K cm -3 . 



In contrast, if the emission enhancement detected by XMM comes from a uniform density 
region with the above properties, the electron density is about 4xl0 -3 cm -3 , which would 
imply a pressure of P/k = 2xl0 4 K cm" 3 . This pressure is nearly two orders of magnitude 
larger than that of the ambient medium. 

A likely mechanism for the production of this emission is the differential motion between 
the HI cloud and the hot ambient medium. This differential velocity is about 300 km s" 1 
( Ivan der Marel et al.l 12002). although recent obser vations show that it may be closer to 370 
km s _1 ( IKallivayalil et al.ll2006l ; iBesla et al.l 120071 ). To be conservative, we will assume that 
the hot ambient medium is at rest with respect to the Milky Way and that the Magellanic 
Stream is moving through it at 300 km s" 1 . As a crude measure of the energy available, we 
estimate the ram pressure heating rate by the expression \pr 2 c v 3 , where r c is the cloud radius, 
p is the ambient gas density and v is the relative velocity. Using this, we find the power 



available is lxl0 35 (n/10" 



cm 



5 ) erg s 1 . This is about a factor of twenty larger than the 



observed power, so there is an adequate power supply to account for the observed radiative 
emission, even if the density is 2.5xl0~ 5 cm -3 . We now consider how this energy might be 
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transferred from relative motion to thermal energy. 

This cloud motion might drive a weak shock into the ambient medium, since the dif- 
ferential velocity is comparable to the sound speed. Even if it were a strong shock, it could 
only raise the density by a factor of four, rather than the factor of about 40 that is observed. 
A more likely explanation is that the emission occurs from the shock that is driven into the 
neutral cloud, heating the gas and causing this overpressure medium to flow away from the 
HI cloud. Additional proce sses are likely to be important, such as mixing between the HI 



cloud and the hot medium (B land-Hawthorn et al. 



20071 : IParriott fc Bregmanl 120081 ) . 



These energy budget and plausibility arguments are helpful in understanding the pro- 
cesses that lead to X-ray emission, but they are no substit ute for a serious simula t ion. S uch 



a numerical hydrodynamic simulation was carried out by iBland-Hawthorn et al.l (120071 ) for 
the purpose of modeling the optical Ha from the Magellanic Stream clouds. In the future, 
it may be possible to extend these calculations to understand whether the interaction be- 
tween a Magellanic Stream cloud and a hot ambient medium can cause the observed X-ray 
brightening. 

We are especially grateful to Robert Braun and Butler Burton for making the H I 
total intensity map of CHVC 125+21-207 available to us. Also, we would like to thank 
Renato Dupke, Chris Mullis, and Steve Snowden for their suggestions and assistance and to 
the Chandra and XMM-Newton teams for answering many of our questions. We gratefully 
acknowledge financial support for this research, which was provided by NASA. 
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Table 1. X-Ray Observations 



Target 


Instrument 


ObsID 


RA (J2000) 


DEC (J2000) 


^>exp (kscc) 


^clean (kscc) 


NGC 891 


Chandra 


794 


02 22 33.4 


42 20 57.0 


51.56 


36.15 


NGC 891 


Chandra 


4613 


02 22 31.3 


42 19 57.3 


120.4 


108.5 


CHVC125 


Chandra 


2253 


12 29 10.0 


75 23 45.0 


48.01 


45.27 


CHVC125 


Chandra 


2484 


12 28 20.0 


75 23 45.0 


24.34 


10.09 


CHVC125 


Chandra 


2486 


12 27 30.0 


75 23 45.0 


22.96 


21.2 


MS30.7 


Chandra 


5038 


00 13 18.9 


-27 13 24.0 


50.14 


49.51 


MS30.7.off 


Chandra 


5039 


00 15 05.5 


-27 27 34.7 


49.88 


33.49 


MS30.7 


XMM 


204670101 


00 12 56.3 


-27 12 06.7 


50.03 


37.0 


NGC 5907 


XMM 


145190101 


15 15 59.0 


56 20 34.4 


52.19 


5.9 


NGC 5907 


XMM 


145190201 


15 15 58.9 


56 20 33.7 


52.74 


6.6 
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Table 2. NGC5907 Measurements 



Method 


Region 


On-galaxy 


Off-galaxy 


On-Off 






(cts ksec -1 arcmkr 2 ) 


(cts ksec -1 arcmin -2 ) 


(cts ksec -1 arcmin -2 ) 


x 2 


1 


1.12 ±0.13 


1.57 ±0.16 


-0.45 ±0.21 


x 2 


2 


1.78 ±0.17 


1.49 ±0.16 


±0.29 ±0.23 


CFW 


1 


1.56 ±0.23 


2.15 ±0.26 


-0.60 ±0.34 


CFW 


2 


2.36 ±0.28 


1.99 ±0.26 


±0.34 ±0.39 



Table 3. MS30.7 Measurements 



Region Flux On-Off Flux 

(cts ksec -1 arcmin~ 2 ) (cts ksec -1 arcmin -2 ) 



1 


2 


.39 ±0. 


15 


0.69 ±0.15 


2 


2 


.66 ±0. 


.15 


0.96 ±0.15 


3 


2 


.56 ±0. 


.15 


0.85 ±0.15 


4 


2 


.15 ±0. 


.14 


0.45 ±0.15 


5 


2 


.27 ±0. 


.14 


0.56 ±0.15 


6 


1 


.71 ±0. 


.12 


0.00 ±0.13 


7 


1 


.95 ±0. 


13 


0.25 ±0.14 


8 


2 


.32 ±0. 


14 


0.61 ±0.15 


9 


1 


.79 ±0. 


.13 


0.09 ±0.14 


10 


1 


.76 ±0. 


12 


0.05 ±0.14 


11 


1 


.84 ±0. 


13 


0.13 ±0.15 


12 


2 


.28 ±0. 


.15 


0.57 ±0.16 


13 


2 


.10 ±0. 


13 


0.39 ±0.14 


14 


2 


.09 ±0. 


13 


0.38 ±0.14 


15 


1 


.90 ±0. 


.12 


0.19 ±0.14 


16 


1 


.85 ±0. 


.12 


0.14 ±0.14 


17 


1 


.92 ±0. 


13 


0.21 ±0.14 


18 


1 


.92 ±0. 


13 


0.21 ±0.15 


19 


1 


.89 ±0. 


13 


0.18 ±0.14 


20 


1 


.96 ±0. 


.14 


0.25 ±0.15 



Table 4. Surface Brightness Contrasts 



Object 


Observatory 


On Cloud 


Off Cloud 


Instrumental 


Unresolved 


On 


Off 


Ratiol 


Local 


Ratio2 






(raw) 


(raw) 


Background 


Pt. Src. 


(net) 


(net) 




Bubble 




MS30.7 


Chandra 


1.29 ±0.042 


1.20 ±0.04 


0.61 


0.04 


0.64 


0.55 


1.16 ±0.097 


0.15 


1.22 ±0.130 


CHVC 125 


Chandra 


1.20 ±0.033 


1.14 ±0.031 


0.61 


0.03 


0.56 


0.50 


1.12 ±0.085 


0.13 


1.16 ±0.114 


MS30.7 


XMM 


3.24 ±0.09 


2.60 ±0.05 


0.70 


0.21 


2.33 


1.69 


1.38 ±0.051 


0.50 


1.54 ±0.068 


NGC 5907 


XMM 


1.90 ±0.11 


1.98 ±0.11 


0.45 


0.25 


1.20 


1.28 


0.94 ±0.125 


0.31 


0.92 ±0.167 


MS30.7 


Sum 














1.33 ±0.045 




1.47 ±0.060 



Note. — The units of surface brightness are cts ksec -1 arcmin -2 , which apply to columns 3-8, and 10. The surface brightness due to the instrumental 
background (column 5) and unresolved point sources (column 6) are removed from the raw surface brightness to produce the net counts (columns 7, 8). 
The quantity Ratiol is the ratio of the net counts on divided by the net counts off the cloud (column 7 divided by column 8). The quantity Ratio2 
(column 11) is the same ratio, after the removal of the surface brightness contribution from the Local Hot Bubble (column 10). The sum of the MS30.7 
results is a weighted average. 
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Fig. 1. — The 21 cm total intensity map of the LMC, SMC, and the Magellanic Stream 
(IPutman et al.l 120031 ) . with the location of cloud MS30. 7-81. 4-118 marked with a circle of 
radius 2°. The intensity values are on a logarithmic scale with black corresponding to 
Nhi > 6xl0 20 cm~ 2 , and the faintest levels corresponding to ~ 2xl0 18 cmr 2 . 
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Fig. 2.— The 21 cm total intensity map of MS30. 7-81.4-118 with the Chandra ACIS-S3 field 
of view (8' square) and the approximate XMM-Newton EPIC field of view (larger polygon) 
superimposed. The four polygons within the ACIS-S3 field show the high and low column 
density regions used for the extraction of the X-ray counts. The scale at the bottom is in 
units of 10 19 cm -2 . In the region of the X-ray observations, the peak H I column density 
is 4xl0 20 cm" 2 , about three times higher than the foreground Galactic column density of 
1.5xl0 20 cm" 2 . 
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Fig. 3.— The 21 cm total intensity map of CHVC 125+41-207 with the Chandra ACIS-S3 
field of view (8' square) superimposed for ObsID 2253, 2484, and 2486 (left to right); the 

In this interferometric H I column density map 



scale at the bottom is in units of 10 18 cm 2 



,20 cm -2 



( iBraun fc Burton! l2000l ). the peak cloud intensity within the ACIS-S3 field is 5xlO z 
while the regions adjacent to the cloud have a mean value of about 3xl0 20 cm -2 . The circle 
marks the location and approximate size of the diffuse enhancement in the X-ray emission. 
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contour lines: 1. 2. 4. 10. 20. 40. 100 * 10 1 " cm 




12 h 30 m 12 h 20 m 12 h 10 m 12 h m ll h 50 

Right Ascension (J2000) 



Fig. 4. — The single-dish 21 cm tot al intensity map o f CHVC 125+41—207 shows that the 
cloud has a head-tail configuration (IBriins et al.ll200ll ). The highest column density regions 
are observed in the X-ray band (see Figure 3). The location of the diffuse X-ray enhancement 
is shown with a white X (also shown above). 
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Fig. 5. — Two Chandra images of NGC 891 in the energy band 0.4-1.0 keV, where the region 
of high optical extinction is delineated by the polygon. Each pixel is 0.5" and the field is 
3.3' high. Each dot is a single photon. The X-ray shadow seen in the first observation (left, 



Bregman fc Irwin! 120021 ) is not confirmed with more recent data (right). 
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Fig. 6. — The H I contours are superimpo sed on the optical image of the edge-on galaxy 
NGC 5907 (jSancisil Il976l : IShang et al.lll998l ). The region where the H I extends beyond the 
optical galaxy is where we have sought to search for a shadow in the diffuse X-ray background 
that is produced from absorption by the cold gas. 



- 29 - 





4 




3 


If) 








c 

D 


2 


O 




O 





200 



300 



400 



500 



600 



700 



800 



Pixels 



0.02 



w 0.00 

c 

o -0.02 



-0.04 




200 300 400 500 600 

Pixels 



700 



800 



Fig. 7. — Horizontal cross-section of the flattened 0.4-1.0 keV image of NGC5907. Both 
panels show row 510 near the center of the image; the top panel covers the full range of 
the positive pixel values, wherease the bottom panel focusses on the pixels that did not 
receive any photon events and now posess negative values after the flattening. The dotted 
line represents zero counts for easier reference. Bad pixels are masked and set to zero after 
the flattening. 
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Fig. 8.— XMM-PN 0.4-1.0 keV image of NGC5907. A bad pixel mask has been applied 
setting correspond ing; pixel values to zero. The shown H I contours (0.69 and 75 mJy beam -1 , 



Shang et al.lll998l ) represent the optical disk of the galaxy and the H I halo. The two pairs 



of rectangular regions used for the X-ray background measurements are shown as well. 



Fig. 9.— XMM-PN 0.4-1.0 keV image of the Magellanic Stream Cloud MS30.7-81.4-118. A 
bad pixel mask has been applied setting correspon ding; pixel values to zero. Only one H I 
contour (iV(HI) = 3 x 10 20 cm~ 2 . lBruns et al.ll2005l ) is shown. The H I column density drops 



sharply on the southwest side of the cloud, but decreases only slowly toward the northeast. 
The twenty measurement regions are labeled and follow the exposure map contours. 
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Fig. 10.— XMM-PN measurements of the Magellanic Stream Cloud MS30. 7-81.4-118. The 
left panel shows the X-ray fluxes in cts ksec -1 arcmin -2 relative to region 6 (whose flux is 
thus zero). The right panel shows the resulting signal-to- noise. The black outline in either 
panel represents the H I contour iV(HI) = 3 x 10 20 cm -2 (IBriins et al.l 120051 ) . 



